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Abstract 


A finite difference computer model is developed to simulate the exposure statistics of a radio 
frequency buoyait aitenna as it is towed in a random seaway. The model allows the user 
to prescribe antenna properties (length, diameter, density, etc.), sea conditions (significant 
wave height, development of sea), and tow speed. The model then simulates the antenna-sea 
luteraction for the desired duration to collect statistics relating to antenna performance. The 
model provides desigii engineers with a tool to predict antenna performance trends, and conduct 
desigu tradeoff studies. The antenna envisioned is a submarine floating antenna which would 
enable communications at speed and depth, greatly enhancing the stealth and survivability of 
the US Navy’s submarine force. 
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Chapter 1 


Introduction 


Driven largely by demand in the commercial sector. the capability for high speed data transfer 
has grown at an amazing rate over the last few years. Recently deploved military coimiuni- 
cations systems have been able to leverage off their commercial counterparts to some extent, 
and now operate at higher frequencies and data rates than the systems they have replaced. 
However, the United States subinarine force faces unique challenges in fully utilizing systenis 
which operate at higher frequencies. Submarine operations are typically driven by the need 
for stealth, and stealth is best preserved by remaining fully submerged. Currently submarine 
comuunications can be typified by a pattern of coming to periscope depth periodically to raise 
a radio mast, send or receive required message traffic. then submerge. While submerged. radio 
comuiunications systems are limited to low data rate receive-only modes transmitted at very 
low frequencies (ELF or VLF bands). Clearly, the communications goal of the submarine force 
is to be able to utilize the high frequency (UHF, SHF, or EHF bands) and high data rate 
colmunication circuits while retaining the stealth of operation at depth. 

This goal of high data rate communication at depth was the motivation of a U. S. Navy Ad- 
vanced Technology Demonstration (ATD). Initially funded jointly by DARPA and ONR in FY 
96-97 as a feasibility study, the Buoyant Cable Array Antenna research became a DARPA/ONR 
Demonstration in FY 98. The program was promoted to the level of Navy ATD in FY 00. This 
ATD. one of only 6 ATD’s funded by the Navy for fiscal year 2000, reflects the high prioritv 


the Navy places on submarine connectivity. The goals of the program include |1] 
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e verify the capability to conduct beamforining signal processing 
e verify adequate antenna performance and gail 11 a seaway 


develop transmit capability 


conduct design studies on key antenna characteristics (diameter, length. density. etc.) 


In addition to the engineering challenges inherent in any new design work, the program 


llanager must ensure the desiguers keep system affordability and flexibility im mind. 


1.1 Background 


Lincoln Laboratory is a federally funded laboratory administered by Massachusetts Institute 
of Technology. Lincoln Lab, specializing i radar and electrical engineering. 1s well qualified to 
tackle the challenge of designing a floating linear array RF communication antenna. However, to 
better characterize the behavior of the antenna in a random seaway, Lincoln Lab added Professor 
Jerome Milgram of the MIT Ocean Engineering department to the antenna development effort. 
Professor Milgram was made the head of the Exposure Statistics Working Group, whose charter 


was to develop a model to provide statistics on 


e the average fraction of antenna elements exposed over an extended period of time 


e the temporal variation of the exposure statistics, cluding duration of reception dropout 


t111es 


e the spatial distribution of those exposed elements 


This thesis work was undertaken to support the charter of the Exposure Statistics Working 
Group. By constructing a finite differeuce computer model of the autenna, a desigu tool was 
envisioned that would allow engineers to evaluate perforniance trends and pose “what if” desigu 
questions. Aid most unportantly, the computer inodel would permit performance trade-off 


studies to be conducted at the system engineering level. 
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Chapter 2 


Model Development 


2.1 Characterization of Sea 


2.1.1 Linear Theory Ocean Waves 


For the purposes of this work. the sea was assumed to be described by linear theory. random. 
deep-water waves. “Deep water” is defined to exist when A < 2d where d 1s the water depth. 


The random aspect of the waves implies that 


e the distribution of wave elevations from the calm water level follow a gaussian forin, and 


have zero 1ea11 


e the sea statistics obtained from a single location are the same as the statistics obtaimed 


from an ensemble of locations (ergodic) 


e the sea statistics remain constant with respect to time (steady state process). 


Because the last condition is so restrictive, we relax it to allow the statistics to vary slowly 
with time (weakly steady state process). A summary of the deep water linear theory wave 
relatious is provided in Table 2.1. 

Assunling a lower sea spectral energy bound of 0.7 rad/sec, which is reasonable for the 
rauge of seas analyzed im this report, this translates into a deep water depth of 252 meters. 
For depths shallower than about 250 ineters, the deep water theory begins to break down, aid 


a depth-dependent treatment of the sea waves must be considered. Referring to the formulas 
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Quantity Formula 
i a a oy °° ©. ——- 

Wave Number — Fi 

1 _— 4a 
Circular Frequency a 
Dispersion Relation — 7A 
Velocity Potential o= A eke sin(ka’ — wt) 
Surface Displacement h = acos(k2z’ — wt) 
Radian Frequency ud = gk 
Phase Velocity Gay oh 
Group Velocity f= gk 
Horizontal Fluid Velocity u = awe" cos(kx’ — wt) 
Vertical Fluid Velocity w = awe™ sin(ka’ — wt) 
Wave Energy Density —— 5 pga° 


Note: 2’ refers to a stationary x coordinate systei. 


Table 2.1: Linear Theory Deep Water Wave Relations 


in Table 2.1, the wave forins in linear theory can be seen to be sinusoidal. However. siniple 


observation of the sea surface shows that the wave forms display a complicated structure which 


has random features. ‘This observed behavior can be approximated by adding together a series of 


slusolds. each frequency component having the appropriate amplitude and a randoi phase. As 


the number of frequency components increases, the approximation approaches the contimuous 


distribution of the actual sea. This discrete method of describing a one dimensional random 


seaway 1s presented in equation 2.1 


N 
ect) = wy a,(w) cos(k;(w)a” — w;t + ¢;) (245) 
t=] 


where 


h(x’,t) is the surface elevation of the sea, referenced to the calm water level 


a;(w) is the wave amplitude of each discrete frequency component 


k,(w) is the wave nuinber of each frequency component 


w, 18 the set of discrete circular frequencies chosen 


@, 1s the random phase assigned to each frequency component 


and the e** term has been dropped, since evaluation occurs at the sea surface c= 0 


It is fairly straightforward to deteriiine most of the terms on the right haid side of equation 


2.1: A; is deteriined from the dispersion relation, ¢; is a random number with a uniform 
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Figure 2-1: Sample Ocean Wave Energy Density Spectrum: H, = 2.0 m 


distribution [0,27]. The evaluation of the ainplitudes, however, requires some discussion. 

We begin by defining the spectral energy density function, S(w), which represents the energy 
of the ocean waves having frequency w. (See Figure 2-1 for a typical energy density spectrum. ) 
The units of S(w) are defined to be in m*/sec, so that when S(w) is integrated with respect 
to w, the result is energy in units of m*. (All that is lacking to give this quantity the normal 
units of energy is the 5g term seen in Table 2.1 for the wave energy density.) The spectrum of 
Figure 2-1 is said to be single-sided, which refers to the fact that the frequeicies are constrained 
to be positive. For a double-sided spectrum (positive and negative frequencies allowed, aid the 
spectruin defined such that S(w) = S(—w)) of equal sea severity, one must divide the spectrum 


by 2, since for constant sea severity the area under the spectral energy density curve must be 
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Figure 2-2: Determination of Component Amplitudes a;(w) 
constant. The amplitudes a;(w) are obtained from Eq. 2.2 
ae | 
—a;(w)* = S(w)Aw (2.2) 


2 


where 

a;(w) is the amplitude of the band of sinusoid waves Aw with center frequency w. 

Eq. 2.2 is really nothing more than a statement that the square of the wave amplitude is a 
lneasure of the average energy of the waves. because both sides of the equation represent energy 
in units of m*. A method of obtaining the amplitudes a;(w) is shown in Figure 2-2. 

Although S(w) has been presented in terms of wave energy content at a given frequency. it 
can be shown froin the Weiner-Kintchine theorem that S(w) is related to the wave amplitude 


record at a fixed location in the ocean. Following Ochi, [2] we define a wave record at a given 
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location as x(t). Then the autocorrelation of «’(t) is given by 


Fa — en) = | (t+ 7) dt (273) 


T— 0 


and by the Weiner-Kintchine theorem, the spectral energy density function, S(w), and the 


autocorrelation function, R(T). form a Fourier transform pair. 


(Ce -f- een AT (2.4) 


Rr) = 5 [swede (2.5) 


This method’s value lies in its ability to convert observed wave records into spectral energy 


density functions. 


2.1.2 Spectral Energy Density Functional Forms Assumed for the Model 


Numerous forms have been proposed by researchers attempting to characterize the sea wave 
environment [3]{4][5][6][7]. The model incorporates the spectral formulation as a subroutine. 
hence it is completely arbitrary from the model viewpoint how the spectrum is defined. In order 
to limit the scope of the simulation parameters, yet retain the ability to vary the sea energy 
form, two wave spectral energy subroutines were developed for the model. Though they may 
appear very different at first glance, most spectral formulations share the basic framework of 
BG. 2.6 
S(a) = Fe Blt (2.6) 
WwW) 

where 

A, B are constants to be determined 

H, is the significant wave height. which is defined to be the average wave height of the one 
third highest waves. Sea severity is most commonly given in terms of significant wave height 
along, researchers. However operators (such as U.S. Navy personnel) are frequently more 
fannhar with sea state. (Appendix E provides a conversion between the two systems.) The two 
spectral formulations chosen for the model are the Bretschneider two parameter spectrum. aud 


the Ochi six parameter spectruin, given in Eqs. 2.7 and 2.8 respectively [2]. 
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2 oreo eee 
S(w) = 7 =i ie las (2.7) 
Oyama a NX 2 | 
S(w) = a S- (eS eS (2.8) 
4 oo T(A;) ita; +1 7 


Bretschneider Two Parameter Spectrum 


The two paraiieters of the Bretschneider spectrum (Eq. 2.7) are significant wave height. H,. aud 
modal frequency, w,. Modal frequency is that frequency at which S(w) reaches its inaximum 
value. aud significant wave height is as previously defined. For the special case of a fully 
developed sea (also referred to as a fully arisen sea). i.e. one in which the wind inparts energy 
to the waves at a rate equal to that at which viscous damping dissipates their energy and wave- 
wave interactions transit energy from one frequency to another; a relationship exists between 


Wm and H,. Designating this modal frequency as wm,,. we have 


g 
Wmpsa — 0.4 fe (2.9) 


which establishes the energy equilibrium characteristic of a fully arisen sea. The two require- 


ments for a fully arisen sea are 


e sufficient fetch (i.e. the distance over which the wind is blowing over the water.) 


e sufficient duration (i.e. the length of time the wind has been blowing.) 


If the fetch or duration are insufficient, the modal frequency will be greater than that given 
by Eq. 2.9; if the winds are subsiding, the modal frequency will be less than that given by 
Eq. 2.9. This is best understood by considering the mechanism by which waves are generated. 
The wind causes ripples (capillary waves) on the sea surface, which is a short wavelength. high 
frequency phenomenon. Thus as a storm is building, wave energy input occurs at high frequency, 
and the energy gradually is transferred through wave-wave interactions to lower frequencies. 
Likewise. as a storm is subsiding, the high frequency waves attenuate most rapidly due to 


viscous effects, and the result is a skewing of the spectral energy toward lower frequencies. By 
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Parameter Value 


en 0.844, 

Ts (oie = 

a ier) 70! 
Oo 1.15¢e70-.039Hs 
d| 3.00 

do 1.54e—0:062H; 


Table 2.2: Ochi Six Parameter Spectrum: Most Probable Values in Terms of Hs 


introducing a term we will call the development Dev. we can quantify the degree of sea arousal 
as piven in Eq. 2.10. 


Wry = Dev -wWmp, (2.10) 


Using this convention, values of Dev < 1.0 correspond to a decaying sea. and values of Dev > 
1.0 correspond to a building sea. Figures 2-3 and 2-4 show representative plots of S(w) as 
development aud significant wave height, respectively. are varied. Note that in Fig. 2-3. that 
while changing the development shifts the location of peak spectral energy density. it keeps the 


area under the curve (and thus the total energy) constant. 


Ochi Six Parameter Spectrum 


The Ochi six parameter spectrum is a generalization of the Bretschneider two parameter spec- 


truni in two respects: 


e two spectral peaks are present, each with their own modal frequency and significant wave 


height. and their contributions add together 


e a shape factor. A, can be applied to each spectral peak to skew the curve right or left. 


independent. of the significant wave height or modal frequency. 


In the form given in Eq. 2.8. each parameter is an independent variable. However. by 
aulalyzing 800 wave spectra from the North Atlantic [2]. one can obtain a relationship for the 
SIX paraineters 11 terms of the measured significant wave height. Ochi reports the results of this 
analysis for a fainily of spectra having a 95% confidence, however the results given in Table 2.2 


represent a single “most probable” spectrum [2]. The model implemented the most probable 
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Figure 2-3: Bretschneider S(w) for Varying Degrees of Development 
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Figure 2-5: Ochi Six Paraineter “Most Probable” Spectrum 


form for obtaining the six paraineters of the Ochi spectrum. A sample plot of the Ochi most 


probable spectrum, showing its characteristic bimodal shape, is provided in Figure 2-5. 


2.2 Coordinate System Description 


The coordinate system adopted for the model is a two-dimensional (x-z) system whose origin 
translates in the x direction at the antenna tow speed U. For the case of head seas, U is defined 
to be negative and the waves travel in the positive x direction. Conversely, for the stern sea case 
(7 is positive and the waves again travel in the positive x direction. Figures 2-6 and 2-7 depict 
the coordinate system employed and the basic antenna configuration for both the head and 
sterl sea cases. The equations listed in Table 2.1 apply to a stationary coordinate system. 


The transformation to a moving coordinate system eimploys the simple relatiouship 


m= x tub (2.11) 
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where 
zr’ is the earth-fixed (stationary) coordinate system 
ris the moving coordinate system whose origin translates with velocity U (as shown im 


Figure 2-6 and 2-7). 


2.3 Forces Acting on Antenna 


The vertical forces acting on the antenna are given m the right hand side of Eq. 2.12 


dv Or dw O1¢ 
m— = T(x, fees Cone vile veo enone ar LT d 


-_— ZA 
dt Ox? dt Ox" ( ) 


where 

mis the mass per unit length of the antenna. assumed uniform 

vis antenna vertical velocity (i.e. a) 

T(x.t) is the tension along the length of the antenna, arising from skin frictional drag 
between the subinerged portion of the antenna and the water 

q is the z coordinate of the antenna central axis 

C'p is the normal drag coefficient of a cylindrical antenna section as it moves in the z 
direction through the water 

w is the vertical fluid velocity 

b(.r,t) is the buoyancy force, which is the resultant of its weight (downward force) and water 
displaced (upward force) 

Cara is the coefficient of added mass for a cylindrical section in a flow field normal to its 
long axis 

Pwater 8 the density of water (for seawater of 15° C, pyater = 1025.9 kg/m?) 

ris the antenna radius 

FE'l is the product of the antenna cross sectional modulus of elasticity and area moment of 
inertia (thus a measure of lateral stiffness). 

Thus, the first term on the right hand side of Eq. 2.12 represents a restoring force charac- 
teristic of a string displaced from its (taut) equilibrium: condition. The tension im the antenna 


is calculated from the time varying wetted surface area along the length of the antenna. A 
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tangential drag coefficient of Cy = 0.0035 was used [3]. (Reynold’s numbers varied from 32.940 
to 494.130 for speed/diaimeter combinations of 3 knots/1 inch to 9 knots/5 inches. respectively. ) 
The second term represents the normal drag force a cvlindrical element would experience as it 
is moved through the fluid with its long axis perpendicular to the direction of motion. A value 
of Cp = 1.0 was used (assumes infinite fluid domain). The fourth term represents the added 
ulass effect experienced by a cylindrical object as the fluid around is accelerating normal to the 
long axis of the cylinder. A value of Casa = 2.0 was selected (assumes infinite fluid domain). 
consistent with standard fluid dynamic texts [9]. The final term represents the restoring force 
due to the stiffuess of the material, and its distribution with respect to the ueutral axis. At 


®) 


this preliminary stage of the antenna design, a material sinmlar to Tygon™ tubing is planned 
for the outer antenna housing. At a nominal design of 3 inch outer diameter. 0.25 inch wall 
thickness. and modulus of elasticity E = 1.5- 10° Pa, we have a relatively flexible antenna of 
Abou = 500 Pani 

It is recogiized that some of these forces have been approximated, namely in that the values 
quoted above are for Cp and Cyj;, in an infinite fluid domain. The antenna operates in close 
proximity to the free surface. and occasionally portions rise completely above the surface of 
the water. A rigorous treatment of the behavior of Cp and Cys, near the free surface was 
not attempted. Rather, a straightforward linear approximation was employed which varied the 
coefficieuts Cp aud Ca;4 from their infinite fluid domain values when the antenna axis was at 
a depth of 7 (i.e. just submerged), to a value of zero when the antenna axis was at a height of 


r (i.e. just breaking out of the water). This dependence of Cp and Cay, on relative antenna 


depth y is shown in Eqs. 2.13 and 2.14. 


y 1 
Cy) = CDirr ie ene at) (2.13) 
y 1 
Cyuraly) = Ci Ainpinte(—5,, + 5) (2.14) 


2.4 Discretization 


The antenna. of length L, is subdivided first into lengths of the physical radio antenna elements. 


These antenna elements, designed at Lincoln Laboratory. are called “jelly roll” antennas because 
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Figure 2-9: “Jelly Roll” Antenna Element Construction {12] 


thev are a conducting sheet rolled up about 13 times. (See Figure 2-9 for the construction of 
an alteuna element.) The element length is designated Lejem. This first level of subdivision Is 
not onlv logical, but necessary for the exposure statistics calculated by the model. As it turns 
out. this type of antenna element performs well as long as it has some sinall amount (e.g. 1 
iin) of its outer housing diameter exposed along its entire length. But if the elenient becomes 
covered by even a small amount of water at any point along its length, antenna performance 
drops dramatically. Thus it is necessary to monitor numerous points within each element to 
deterinine whether a given element is exposed or not. This leads to the final level of aitenna 
subdivisiou: each element is ineshed into some integer nuinber AJ of sections with length As. 
The subdivision process is shown in Figure 2-10. 

An interesting feature of the program version which calculates the sea elevation and fluid 
velocities by Fast Fourier Transforms (FFTs) is that the frequency meshing (Aw) of the spectral 
energy density function S(w), is now related to the spatial meshing (Aw) of the sea, and thus the 
spatial meshing of the antenna (As). This means that there is an extra constraint to consider 


in selecting the finest antenna mesh size As; or conversely. given a mesh size As, coustraints 
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Figure 2-10: Antenna Subdivsion 


26 





are imposed on the Aw meshing of the spectral energy density function S(w). Some freedom is 
regained if one decouples the length of the anteina and the length of the sea smuulated. This 
decoupling is tmplemented in the model. so that for a typical run the antenna may be 25 m 
long with a mesh size of 0.0625 m, while the length of sea modeled is 512 m long with a mesh 
size of 0.0625 m. The key point is that the mesh length between the antenna (As) and the sea 


(Avr) urust match. 


2.5 Computational Method 


2.5.1 Sea Domain Properties 


The sea domain variables of interest are 


e sea elevation. /t 
e fluid horizontal velocity, u 


e fluid vertical velocity, w. 


Two inethods were employed m arriving at the sea domain properties, each with certain 
advantages aud disadvantages. The first method will be called Sine Wave Superposition (SWS). 
and the second method will be called Fast Fourier Transform (FFT). The source code for both 


of these computer programs can be found in reference [10]. 


Sine Wave Superposition (SWS) Method 


As previously inentioned in Eq. 2.1, the sea elevation can be obtained by a superposition of 
sine waves having random phase and amplitude prescribed by Eq. 2.2. In like maimer the fluid 
velocities can be calculated by the superposition method. The forimulas for h, u and w in an 


earth-fixed reference frame are provided in Eqs. 2.15 - 2.17. 


N 

ee S > ai(w) cos(k;(w) x" — w,;t + ¢;) (Za) 
7=1 
N 

a RCAC w cos(k;(w) x — w;t + @;) (2.16) 
i=1 


2/ 


: 
ip 





N 
w(z’,t) = S_ ai(w) wsin(k;(w) 2’ — wit + @;) (2.17) 
eal 


By invoking the coordinate transformation and the dispersion relation (Eq. 2.11 and Table 2.1. 
respectively), one arrives at the forms appropriate to the coordinate system moving with the 


autenua. These are presented in Eqs. 2.18 - 2.20. 





ee a ec (w) cos | k;(w) a — (wa; — S \t+ , (Zales) 
1=] = 
id ie 

eae) S_ a,(w) w cos |kj(w) x — (w, — ae + Od; (218) 
i=] 
a w*U 

[OT a S_ ai(w) wsin |k,(w) 2 — (w; — wae + o; (2.20) 
7=] ; 


The spectral energy deusity function S(w) used im this formulation was a single sided one. The 


advantages of the SWS method are 
e suuplicity m codmg fluid property subroutines 
e good resolution for S(w) at lower frequencies, where most of the wave energy resides 


e frequency meshing (Aw) is unrelated to the spatial meshing (Ax). allowing both antenna 


element length aud wave energy frequency resolution to be arbitrarily chosen. 


Fast Fourier Transform (FFT) Method 


This method takes advantage of the fact that Eqs. 2.15 - 2.17 are a Fourier series, and thus the 
functions h,u,and w are the inverse Fourier transform of some expression. Since the calculations 
to determine the sea domain properties occur at a fixed pomt m time, the transform variables 
are not w and ¢t, rather they become 6 and x. Table 2.3 summarizes the analogous relationships 
between the time-frequency domains and the space-spatial frequency domains. (The first two 
lies may be taken as definitions; the last line utilizes the dispersion relation.) Because the FFT 
algorithm used for the model assuined the presence of both positive and negative frequency 
components, If was necessary to use a two sided spectrum for this method. In order to retain 


the same amount of energy (i.e. area under the spectral energy density curve) as the transition 
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Table 2.3: Comparison of Fourier Transform Domain Relationships 


from one-sided to two-sided spectra is made, it is necessary that the formula for the coimponent 
aplitudes be adjusted from that given m Eq. 2.2. In addition, it 1s convenient to mcorporate 
the random phase 9, into a complex component aniphtude c;. With these changes im mind. the 


formulas providing the sea domain properties are given in Eqs. 2.21 - 2.25. 


h(x) 


Re {oath c;(b;) exp |2 i (27b;2 + (27b,U — V/27gb;) t)| (oF 


u(x) Re {So tt GAO v 2ngbi exp |2 i (27b;x + (27b,U — V/27Qgb;) t ) (2.22) 


w(x) = Re p2 ~—ic;(b;) /27gb; exp |i i (Qrb;x2 + (27b;U — \/27gb,) 9]| (2223) 


| 
Re {eb} = 5 S(bi) Ad cos(o,) (2.24) 
lia One — = S(b ;) Ab sin(¢;) (2.25) 


Then the Fourier transform pairs become 


Jal (2) Ss Ge) (2.26) 
JV2ngb; H(b) & u(x) (2°20) 
—i\/27ngb; H(b) = w(x) (2.23) 


where 
N 


H(b) = S~c,(b) exp it (276 — V2rgb;) (2.29) 


—a 
and 


LD ee /3 (22) (2.30) 
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is defined to meal 


N 
Ath). = 3 h(x) exp (—i27b2) eoh) 


=!) 


H(b) exp (22702) (2332) 
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By setting the complex component amplitudes c; to obey the relation 





Ci(b,) = e(—);) (2.33) 


that is. the positive frequency coniponents are the complex conjugates of the negative frequency 
components. then the niaginary parts of Eqs. 2.21-2.23 are zero, and the Re{ } can be dropped 
from the equations. When converting between expressing the spectral energy density fulction 
as a function of w to another variable, such as b, the basic relations of Eqs. 2.34 aid 2.35 are 


useful. 


S(w)dw = S(b) db (2.34) 
S(b) = Sw) (2.35) 


An important feature of the FFT method is the linkage of the spatial frequency domain 
meshing (Ab), and the sea spatial meshing (Az). As previously mentioned. the key to the FFT 
method is to make the antenna spatial meshing, As, match the sea domain spatial meshing, 
Az. This is accomplished by using the relations given in Eqs. 2.36 - 2.39, and ensuring that 


the sea Az equals the desired antenna As. 








1 | 
Ab = 7 — = bmn (2.36) 
Ie ANG iNeed ] | 
Ce eee eee Dae 
mes 2 QLsea  2Ar Ca 


Wmin = Vy 20 GUmin (2.38) 
Wena =" We OU max (2:30) 


The primary advantage of the FFT method is computational speed. The FFT method 
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Frequency, w; 


Component SWS Method FFT Method 
] 0.347 0.347 

2 0.381 0.491 

5 0.414 0.601 

4 0.448 0.694 

5 0.482 0.776 

6 0.515 0.850 

7 0.549 0.918 

8 0.583 0.981 

9 0.617 | 1.041 

10 0.650 1.097 
Notes: 

1. FFT case: w € (0.347, 22.2] with 4096 frequency 
components 

2. SWS case: w € [0.347, 24.0] with 700 frequency 
compouents 


Table 2.4: Comparison of Frequency Resolution for the SWS and FFT Methods 


executes approxunately 10 tunes faster than the SWS version. A further increase in execution 
speed is accomplished by not computing the FFT every time step, At. Rather, the sea data 
( h,u,w) is computed at nAt intervals, and linearly interpolated in between. The mteger n is 
determined by the computer model using the highest sea energy frequency component. Wmax- 
aud the initial estimate of the time step. At. This ensures that the linearly interpolated sea 
evolves in a sinooth and continuous fashion. 

Its disadvantage is that it has poorer low frequency resolution of the sea wave energy than 
its SWS-counterpart. This is demonstrated by a typical run case in Table 2.4 where the lowest 
10 frequencies defined are listed for both methods.As one can readily calculate froin integrating 
Eq. 2.7 with respect to w. for a sea severity of 2.0m significant wave height, 67.7% of the sea 
wave energy lies below a frequency of 1.097 rad/sec. Thus the FFT method suffers in that it 
places inost of the frequency components at high frequencies. where little of the wave energy 
resides. This undesirable result is a consequence of choosing a 6 space meshing with constant 
Ab. which is necessary to compute FFT’s in the b and x domains. (Constant Ab meshing 


produces a non-constant Aw meshing with Aw; = /27qg)nin (vi —Vi- 1).) 


3] 


a 





2.5.2 Integration Method 


The time marching integration method implemented is an explicit forward Euler schenic. Such 
schemes are characterized by using information at time ¢t to calculate functions at time ¢t — Af. 
A simple example of how the forward Euler method is used is given below. 


Solve 


du 9 

—S> = ), 
7 u (2.40) 
u(0) = 1 (2.41) 


The tine derivative is approximated using the Euler method as 








du pnt — Un 
a 2 (2212) 
dt ae fn+l _ yn 
which for constant time step At becomes 
du prt Saye 
a = (2.43) 
oh ae At 
having adopted the notation of Celia and Gray [11], where 
UO at’) (2.44a) 


aud t” is the n“” time step. Then the solution proceeds step by step, beginning with the initial 


value of u given in Eq. 2.41 
du 


7 = u(0)? = 1 (2.45) 


‘=0 





then solving Eq. 2.43 for U"*t! we obtain 





d 
Uy) = ys At — (2.46) 
dt |4=0 
With the value of u now known at time t (U'), we can calculate 
du 2 
as = (U etn 
dt t=] ( \ 0 
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and 


du | 
ve=Ul sat — (2.48) 
dt t—ti 


In like mammer, we move forward in time to calculate any desired U”. 

Spatial derivatives ( 5-4) alternate between forward centered and backward centered as the 
derivative order n increases. accomplishing a space centered scheme. Thus the algorithin can 
be classified as a forward time, centered space (FTCS) method. For many desigus cousidered. 
the antenna was approximated as a flexible string of negligible stiffMess. With this in mind. Eq. 
2.12 becomes a string vibration type equation with damping. As such the Courant condition 


[11. p. 231] becomes a useful starting pomt in selecting a time step size. The Courant condition. . 


| eres, 
ct 
eee 
co a 


] (2.49) 
where 

c is the speed at which waves (i.e., disturbances) propagate down the antenna 

At is the time step size 


Av is the antenna spatial mesh size 


is used to calculate an initial time increment based on the string wave velocity given by Eq. 


0. 


to 
Cn 


(2.50) 


= 
p 
hiPea. 2:50 

T is the string tension, and 

p is the linear string density. 

Because the tension varies with time, and the niaxiinum tension varies with sea severity 
(due to tangential drag associated with the square of local fluid velocities), the calm water 
tension is used in Eq. 2.50 to form an initial estimate of the time step. Then by multiplying the 
calin water time step by a safety factor, a time step appropriate to the sea severity is arrived 
at. (For a 3.0 1 significant wave height, this safety factor is typically set to 0.25) When the 
autenna 1s modeled as having stiffness. the velocity at which disturbances propagate down the 
autelna is now frequency dependent, and the above iuplementation of the Couraut condition 


no longer applies. Numerical experimentation was used to determine the safety factor which 
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produced au acceptable time step size. 


2.5.3. Boundary Conditions 


The boudary conditions imposed on the anteuna endpoints are given in Eqs. 


q(0) = h(O) (2.51) 
oe =) he) | (2.52) 


That is. the anteiua follows the sea surface at its e1dpoints. This is envisioned to be accom- 
plished at the trailing end by a drogue having positive buoyancy. aud at the tow point by one 


of several options. These options include 
e a lifting body 
e a surface following buoy 
e a leader of positive buoyancy 


e some conibination of the above features. 


The problem pertaining to the antenna leader (that portion from the submarine tow point up 
to the surface) is being studied by another graduate student at MIT. and will not be addressed 
in this thesis. 

2.5.4 Initial Conditions 


The antenna simulation begins at t = 0 with 








QU ea) (2255) 
dq(xz) _— dh(z) : 
aA = Pr (04) 


that is. the antenna is placed at the sea surface, and given ai initial vertical velocity matching 


that of the sea surface. 


34 





Chapter 3 


Verification of Model Seas 


3.1 Significant Wave Height 


The random seaway was constructed by the SWS and FFT methods. as described in Sections 
2.5.land 2.5.1. By running the sea generation subroutines, and sampling the wave amplitude 
record, one can determine the significant wave height of the sample. This can then be conipared 
to the significant wave height which was originally specified. From Ochi[2) we have that the 


sample’s siguificant wave height can be found by Eqs. 3.1 and 3.2, 





il N 
— Soh eat 
- ie (3.1) 
H, = 401-4 (3.2) 


where 

fh; is the sanipled wave amplitude, measured froin the calin water level 

a is the standard deviation of the sampled wave amplitudes 

H, is the measured significant wave height of the sample. 

(The “ ~~” symbol denotes that the quantity is computed from a sainple, and therefore nay 
deviate somewhat from the theoretical value.) 

A sea severity of 3.0 in sigmficant wave height was simulated by both the SWS and the 
FFT methods. Twelve wave ainplitude records were analyzed by sampling the seas at 8.0 m 


intervals. The results are summarized in Table 3.1. Good agreement to the specified significant 





FFT SWS 
Sainple Mean(/; ) a H i Mean(h; ) a H. 
(11) (1m) (in) (in) (m) (i) 
l -.0008 0.7635 30214 -0.0005 0.7463 2.992s 
2 0.0002 0.7512 3.0124 0.0002 0.7498 3.0068 
3 -0.0003 0.7490 3.0037 -0.0006 0.7464 2.9932 
4 -0.0003 0.7481 2.9997 -0.0006 0.7494 3.0051 
5 0.0006 0.7511 3.0119 -0.0010 0.7475 2.9976 
6 -OOTLOmaUsrolon c. 0141) 0.0010 0.7484 3.0009 
7 0.0001 0.7504 3.0092 0.0000 0.7503 3.0088 
8 -0.0001 0; 7e0S es 0112 -0.0002 0.7501 3.0081 
g -0.0006 0.7490 3.0033 0.0018 0.7483 3.0007 
10 -0.0004 0.7489 3.0032 0.0011 Q27525° 3:0l7y 
al 0.0002 0.7492 3.0044 0.0002 0.7479 2.9991 
12 0.0002 0.7500 3.0075 -0.0005 0.7505 3.0095 
Std. Dev.’s: 0.0005 0.0015 0.0061 0.0008 0.0018 0.0072 
Averages: 0.0000 0.7502 3.0085 -0.0001 0.7489 3.0033 


Table 3.1: Simulated Significant Wave Heights for a Specified Sea Severity of 3.0 m 


wave height is achieved. with an average overshoot in H, of 0.28% and 0.11% for the FFT and 


SWS inethods. respectively. 


3.2 Wave Period 


Determination of the wave period (or wavelength, frequency. etc.) becomes difficult for a 
randoin seaway, where wave shapes and heights are constantly evolving. However, it is quite 
uecessary to verify the wave period, because it is possible to correctly model the sea’s significant 
wave height, but incorrectly model the wave period (which the research inadvertently proved). 
The inethod used to verify the wave period iivoked the superposition principle: the sea waves 
generated by both methods are ultimately a superposition of individual sine waves, thus the 
verification of wave period was restricted to seas represented by a single harmonic wave. An 
Muportaut aspect of this process is to verify that the antenna sees the correct “encouuter fre- 
au effective change of the wave frequency as the antenna is towed in a seaway. 


quency, 1.e. 


The encounter frequency, we, is given by 


2 
Wy is OG 
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Figure 3-1: Variation of Encounter Frequency vs. Tow Speed 
where 


wo is the wave frequency in a stationary coordinate system 

U is the antenna tow speed 

g is the acceleration due to gravity. 

Thus for towing into head seas ( U is negative), the encounter frequency goes up, and for 
towing with stern seas ( U is positive), the encounter frequency goes down. The encounter 
frequency was measured for a harmonic wave of frequency wo = 1.0 rad/sec at tow speeds from 
-6.0 un/s (towing into head seas) to 6.0 m/s (towing with stern seas). The resulting plot of we 
vs. U is shown in Figure 3-1. The plot of we vs. U is expected to be linear with a slope of = 
(see Eq. 3.3), which for wo = 1.0 rad/sec has a value of -0.1019 m7!. The measured slope was 
-0.1017 m~!, which shows the proper dependence of encounter frequency on speed has been 


captured. 





Chapter 4 


Exposure Statistics 


The primary goal of the effort was to calculate exposure statistics as various design parameters 
were varied. However, before proceeding on to the exposure statistics results. the issue of 
“what constitutes an exposed (i.e., receiving) antenna element?” must be addressed. The 
experimentally observed performance characteristics of the “jelly roll” antenna elements drove 
this definition. It was found that for a level element in calm water the antenna gain is stable 
right up to the point at which the outer skin becomes completely submerged. It was also 
found that if any portion of the outer skin was covered by water (such as might happen in a 
wave field), then the gain of the element suffered drastically. Thus an exposed element was 
defined to be one which was exposed by some small threshold height (set to 0.04 Deameter for 
all simulations) at every point along the length of the element. In a discrete computer model 
some approximation is necessary with regard to “every point along the length of the element.” 
Thus each antenna element was oversampled by an integer factor, AY, which the user prescribes. 


Typical designs considered had an element length of 0.5 m, and were oversampled at A/ = 8. 


4.1 Selection of Simulation Series 


In order to bound the problem, it was necessary to discuss with the customer (U.S. Navy) what 
the operational profile required, and with designers (Lincoln Laboratory) what the engineering 
aspects of the antenna gain required. These discussions led to the following conclusions on the 


design space to be simulated: 
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e tow speed: U < 10 knots 

e sea severity: H, < 3.0 m (sea state 4 or less) 
e length: vary about 2 = 25 m 

e <liaimeter: vary about Diameter = 3 inches 
e density: vary about p = 400 kg/m° 


e antenna element length: vary about Lejem = 0.5 m. 


In addition to perforining simulation series to support the above desigi space, it was neces- 
sary to devote some simulations to verify that the results are, in fact, meaningful and consistent. 


Toward this end snnulation series were made to 


@ collpare exposure statistics results of the FFT and SWS methods 


e vary duration of the simulation to determine what constitutes an adequate length of time. 


It is not clanned that these sunulation series represent an exhaustive list of all the parameters 
which could be varied, nor is it claimed that the model is completely verified by the cousisteucy 
runs made. However the trends presented will be valuable for the design engineers as thev 
consider performance trade-offs, and will allow them to better explore the design space of 


interest. 


4.2 Results 


4.2.1 Statistics Collected by the Model 


The coinputer model was prograinmed to track three key statistics: 


e average exposure — the average fraction of aitenmna elements exposed throughout the entire 


duration of the simulation 


e threshold inet — the average fraction of the time that some specified number of antenna 


elements (or more) were exposed 
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Table 4.1: Exposure Data Matrix File Format 


e average counterdetection length — the average length of the autenna completely exposed 


(i.e.. completely out of the water). 


(The counterdetection length, while not in the origmal list of desired statistics. was added 
because of the tactical significance of compromising submariie stealth.) 

In addition, a record is made at each sample time of the condition of every antenna element. 
such that when the simulation is finished, an exposure matrix is available for further post- 
processing. The structure of the exposure matrix file (an ASCII text file) is shown in Table 4.1. 
(A “1° mdicates an exposed antenna element, a “O” indicates a non-exposed element. The x; 
represent antenna element locations, the t; represent sample times.) Examples of information 


which could be extracted from the exposure matrix include 


e average drop-out time of a specific element 


e average (spatial) length of wave wash-overs 


histograms and standard deviations of the number of exposed elements 


exposure statistics for an array not fully populated (i.e.. one could examine only the 


columus x; which had antenna elements installed) . 
Three output files are generated for each run of the computer model. They are the 


e End-state sunnuary file — lists input parameters used for the sinulation, forces and veloc- 


ities for the final time, aud a summary of statistics collected 
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Parameter Description 

Method SWS or FFT used to generate seas 

jou significant wave height 

Devel sea development factor 

Ga lowest wave frequency component 

Waa highest wave frequericy component 

N freq number of (positive) wave frequency components 
L length of antenna 

Nera number of antenna elements 

es antenna eleinent length 

Al element oversample factor 

As anteimiua inesh size 

Diameter autenna diameter 

Bans arltelina density 

EI anteuna stiffiess 

Duration slinulation duration 

U tow speed 

Nee oe threshold number of elements required for speci- 


fied array gain 


Table 4.2: Simulation Parameter List 


e Exposure data matrix file — provides exposure coudition of each element at each specified 


time increment (refer to Table 4.1) 


e Snapshot file — provides a visual display of the antenna as it floats on the sea surface. 


These output files are described in greater detail 1 Appendix A. 

The input parameters which are specified for each run of the computer model are listed in 
Table 4.2. 

In the-exposure statistics results given below, each series of simulations is presented in the 


following manner: 


e the parameter Input table used in the series is listed 
e plots are provided of the three exposure statistics tracked 
— average fraction of elements exposed 


— fraction of the time that the threshold number of elenients were exposed 


— average counterdetectable length of the antenna 


Al 





— 


ee ee 





It is nnportant to realize that these plots represent averages im time and space. This is 
necessary in order to display the results as functions of key antenna design parameters (length. 
diameter. etc.), which is most useful from a designer’s viewpoint. From the output files. it is also 
possible to analyze how the statistics vary as time progresses. or how the statistics vary down 
the arrav from the tow poimt to the trailing end. Because of the large nuinber of simulation 


runs. only a few representative plots of this type of information will be shown. 
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Parameter Description 
Method SWS 

id various 
Devel 1.0 

ee 0.4 rad/sec 
eae 24.0 rad/sec 
N freq 700 

L 20 Mm 

Netem 30 

ole 0.5 m 

M 8 

eXG 0.0625 m 
Diameter 3.0 inch 

ne 408 kg/m? 
oe 0 Pa in? 
Duration 180 sec 

U various 
Ninresh 40 


Table 4.3: Series 1 Input Parameters 


4.2.2 Series 1: Statistics vs. Sea Severity and Tow Speed 
Input 


The input parameters for Series 1: Statistics vs. Sea Severity and Tow Speed. are given in 
Table 4.3. The tow speeds used were 3. 5, 7, and 9 knots; the significant wave heights used 


seme 0-00 1.0, 1.6, 2.0.2.5. and 3.0 ineters. 


Results 


Average fraction of elements exposed. Figure 4-1 shows a rather weak dependence of 
the average fraction of elements exposed versus significant wave height, decreasing as the sea 
severity increases. However, the plot shows a strong dependence of average exposure versus 
tow speed. 

Fraction of time threshold met. With a required threshold of 40 out of 50 elements. 
Figure 4-2 again shows that tow speed is the dominant factor in the exposure probabilities.. 
The shght rise m threshold met versus significant wave height for the 7 and 9 knot data is an 


uiteresting feature, which will be discussed further below. 
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Average counterdetection length. As one would expect, Figure 4-3 shows the average 
length of the antenna totally exposed increases with tow speed and sea state. The tow speed 
contributes added tension m1 the antenna, and the wave height provides the troughs necessary 


for a section to become exposed. 


Histogram Analysis for Number of Exposed Elements 


As mentioned above, there is slight rise in ‘threshold met’ for 7 and 9 knot runs as the sea 
severity creases. This is not obvious at first, but can be explained as follows. Figures 4-4 - 
4-6 show the histograms of the number of exposed elements over time running at 7 knots in 
significant wave heights of 1.0, 2.0. and 3.0 meters. For low sea states. the waves have little 
energy. and thus have little ability to affect the exposure statistics. This leads to a relatively 
tightly grouped distribution, as in Figure 4-4. Because the mean is below the threshold required. 
the ‘threshold met’ is very low. (As in Figure 4-2.) But as the significant wave height increases. 
the waves contain more energy. and cause greater variation im the exposure statistics. This 
causes the distributions to flatten out. while the mean exposure stays relatively constant. As 
the upper tail of the distribution moves outward from the mean, it provides more “hits” above 


the required threshold of 40 elements. This trend is confirmed in Figures 4-5 and 4-6. 


4.2.3 Series 2: Statistics vs. Length 
Input 


The input parameters for Series 2: Statistics vs. Length, are given in Table 4.4. The threshold 
uumber of elements was varied in proportion to the length (and the nuinber of elements), but 


was kept at 80% of the total number of elements. Element length was kept coustant. 


Results 


Average fraction of elements exposed. The average fraction of elements exposed drops off 
quickly at. first, then becomes nearly flat, as shown in Figure 4-7. This effect must be studied 
more Closely, since the boundary conditions have a much greater relative iupact on a short 
autenna than a long one. (Recall the tow point and trailing end are set to float at the sea 


surface. ) 
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Figure 4-1: Average Exposure vs. Tow Speed and Significant Wave Height 


Parameter 


Method 
Af, 

Devel 

W min 

max 

IN freq 

iE 

ele. 
elem 

Al 

As 
Diameter 
Pant 

eal 
Duration 
U 

pean 


Description 


SWS 

22081 

1G 

0.4 rad/sec 
24.0 rad/sec 
700 

various 
various 

0.61 m 

10 

0.061 in 

3.0 inch 

408 kg/m? 
0 Pa in4 

180 sec 

oO knots 
various (0.8 Netem) 


Table 4.4: Series 2 Input Parameters 
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Figure 4-2: Threshold Met vs. Tow Speed and Sigmificaut Wave Height 
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Figure 4-3: Counterdetection Length vs. Tow Speed and Significant Wave Height 
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Figure 4-4: Histogram of Exposed Number of Elements at 7 kts: H,= 1.0 m 
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Figure 4-5: Histogram of Exposed Number of Elements at 7 kts: H,; = 2.0 m 
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Figure 4-6: Histogram of Number of Exposed Elements at 7 kts: H,= 3.0 in 








Parameter Description 


Method FFT 

ele 27071) 

Devel 1G, 

es 0.347 rad/sec 
ee 22.206 rad/sec 
ee 4096 

iz 20 m 

Netem various 

fe een varlous 

Al varlous 

NS 0.0625 m 
Diameter 3.0 inch 

DES 408 kg/m® 

EI 0 Pa m? 
Duration 180 sec 

U 5 knots 

IN aeeah various (0.8 Netem) 


Table 4.5: Series 3 Input Parameters 


Fraction of time threshold met. Because the threshold was kept at 80% of the total 
number of elements, we see a sharp drop off in the threshold met versus length in Figure 4-8. 
If the threshold were maintained at some constant number of elements (e.g., 15), one would see 
the opposite trend as length increased. 

Average counterdetection length. The rise in counterdetection length in Figure 4-9 is a 
colnbination of two factors: maximuin tension grows proportional to the length, and as antenna 


length grows beyoud the average wavelength it bridges the wave peaks more often. 
4.2.4 Series 3: Statistics vs. Element Length 


Input 


The input parameters for Series 3: Statistics vs. Element Length, are given in Table 4.5. The 
total antenna length was held at 25 meters, and the element lengths used were 0.25, 0.5. 1.0. 
aud 5.0 meters. ‘The oversample factor AJ was varied to keep the antenna mesh length As 


constant at 0.0625 meters. 
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Figure 4-7: Average Exposure vs. Length 
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Figure 4-8: Threshold Met vs. Length 
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Figure 4-9: Counterdetection Length vs. Length 
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Results 


Average fraction of elements exposed. The results m Figure 4-10 show the effect of 
requiring that a contiguous leugth (the element length) be exposed. As the coutiguously exposed 
length requirement grows, the probability of meeting it decreases. 

Fraction of time threshold met. Figure 4-11 shows that the threshold met drops drasti- 
cally as element length grows. The reason can be traced back to Figure 4-10. which shows the 
mean exposure falls below the threshold of 80%. The more the average exposure falls below the 
required threshold number of elements, the ore rapidly the ’threshold met’ drops off toward 
zero. This pattern would be expected for any distribution with central tendencies. (“Ceutral 
tendency distributious are distributions which have a higher probability of returnimg a value 
near the mean than values far away from the mean. The normal (bell-shaped) distribution is a 
classic example of a distribution with a central tendency.) 

Average counterdetection length. Varying the element length should have uo effect on 
the totally exposed length (the elements are assuined to have flexibility given by EJ). Figure 


4-12 supports this claim. 


4.2.5 Series 4: Statistics vs. Diameter 
Input 


The input parameters for Series 4: Statistics vs. Diameter. are given in Table 4.6. All para- 
meters were held constant, except for the antenna outer diameter, which took on values of 1.0. 


20) 30. and 5.0 inches: 


Results 


Average fraction of elements exposed. The average fraction of elements exposed is a 
strong function of anteina diameter. as seen im Figure 4-13. This is to be expected, since 
antenna tension aud normal drag forces (per unit length) grow proportional to the diameter. 
and buoyaucy forces (per unit length) grow as diameter squared. 

Fraction of time threshold met. (Figure ??) Typical for distributions with central 


tendeiicies, and a inean which varies as im Figure 4-13. 
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Figure 4-10: Average Exposure vs. Eleinent Length 
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Figure 4-11: Threshold Met vs. Element Length 
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Figure 4-12: Counterdetection Length vs. Element Length 





Parameter Description 
Method FFT 

ide 2.0t 

Devel eG 

arn 0.347 rad/sec 
Ce 22.206 rad/sec 
ee 4096 

iL Zo 11 

Netem o0 

ae 0.5 m 

Ml 8 

As 0.0625 m 
Diameter various 

Dont 408 kg/m? 
es 0 Pa m4 
Duration 180 sec 

U 5 knots 
Nthresh 40) 


Table 4.6: Series 4 Input Paraineters 


Average counterdetection length. Figure 4-15 shows a strong dependence of counter- 
detection length on diameter. This is again because the buoyancy forces depend on diaimeter 


squared. and thus large diameter desigis “ride” the waves, instead of cutting through them. 
4.2.6 Series 5: Statistics vs. Density 


Input 


The input parameters for Series 5: Statistics vs. Density, are given in Table 4.7. All parameters 
were held constant, except for the antenna density, which took on values of 200. 300, 408, 500, 


600. aud 700 kg/m*. (The density of the prototype antenna element is 408 kg/m°.) 


Results 


Average fraction of elements exposed. The result shown in Figure 4-16 looks very much 
like an inverse of the result from the previous series. This is because the downward weight 
force (per unit length) grows proportional to the density. while other forces modeled remain 


unaffected as the antenna density varies. 
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Figure 4-13: Average Exposure vs. Diameter 


Parameter 


Method 
Af, 

Devel 
min 

UW max 

N freq 

L 

Nelem 
Lelem 

Af 

As 
Diameter 
Pant 

legs 
Duration 
U 


Nt hresh 


Description 
FFT 


peor 

1.0 

0.347 rad/sec 
22.206 rad/sec 
4096 

oon 

o0 

een 

8 

0.0625 m 

3.0 inch 
varlous 

0 Pa m4 

180 sec 

5 knots 

40 


Table 4.7: Series 5 Input Parameters 
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Figure 4-15: Counterdetection Length vs. Diameter 
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Fraction of time threshold met. (Figure 4-17) Typical for distributions with central 
tendencies. and a mean which varies as in Figure 4-16. 

Average counterdetection length. Figure 4-18 shows a knee in the curve at the 408 
kg/im® density. While other factors shift this knee left or right, the density of the prototype 
autenna element happens to be about the lowest value one would desire from a counterdet ection 
staudpoiit at this speed-sea state combination. So from a connectivity viewpoint, lower density 
is desirable: from a stealth viewpoint, lower density is not desirable. Such trade-offs occur 


routinely in subimnarnne operations. 


4.2.7 Series 6: Statistics vs. Duration 


This series and the one following are intended to support the credibility of the above simulation 
results. The purpose of this validation series was to demonstrate that the 180 second duration 
used in the computer runs was long enough to be representative of a random process. There 


were two conflicting goals with regard to the duration: 


e computer run time dictated the shortest duration feasible 


e statistical confidence dictated the longest duration possible. 


Thus three runs were conducted at three different durations (for a total of mine runs), in 
order to understand the variance involved in the statistics presented im the first 5 series of 
snnulations. 

Input 

The input parameters for Series 6: Statistics vs. Duration, are given in Table 4.8. Three runs 
each were inade for duration values of 60, 180. and 720 seconds. 

Results 


Referring to Figures 4-19 - 4-21, each duration has three bars plotted which correspond to the 
three runs made for that duration. The dashed lines represent the best estimate of the actual 


statistic plotted. since it is the average of the three 720 second run results. Above each duration 
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Figure 4-16: Average Exposure vs. Density 
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Figure 4-17: ‘Threshold Met vs. Density 
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Figure 4-18: Counterdetection Length vs. Density 


Parameter Description 
Method FFT 

sale 2.0 1 

Devel 1.0 

copes 0.347 rad/sec 
W max 22.206 rad/sec 
N freq 4096 

L 207) 

Netem 30 

Lic 0.5 m 

MI 8 

As 0.0625 m 
Diameter 3.0 ich 

nee 408 kg /in° 
rk 0 Pa m4 
Duration varlous 

U 5 knots 

ee 40 


Table 4.8: Series 6 Input Parameters 


value is the niean and standard deviation associated with those three simulations. As expected. 
the longer the simulation lasts. the smaller the standard deviation becomes. The goal of the 
colputer model is to show trends which will assist the designers in choosing the optimum 
autennua paraineters. At this stage of the simulation effort, the precise values for elemeut 
exposure cannot be predicted without more extensive experinental validatiou efforts. Thus the 
standard deviation associated with the 180 second duration was judged to be a satisfactory 


compromise between computer run time and statistical confidence. 


4.2.8 Series 7: Statistics vs. Sea Generation Method 
Input 


This series ran identical parameters (within the limits of the meshing conditions) using the 
SWS and the FFT method. in order to ensure the results were consistent between the methods. 
Table 4.9 shows the values used. and indicates which parameters were not latched exactly 
because of frequency meshing issues. The significant wave height was set to 1.0 and 3.0 meters, 


un order to compare the methods over the full range of sea states simulated. 
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Figure 4-19: Average Exposure vs. Duration 
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Figure 4-20: Threshold Met vs. Duration 
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Figure 4-21: Counterdetection Length vs. Duration 
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Parameter 


Method 
As 

Devel 
Wmin 
Wmax 

N freq 

ie 

Netem 
Lelem 

NM. 

As 
Diameter 
Pant 

El 
Duration 


U 
N; hresh 


Description 


FFT / SWS 
Various 

0 

0.347 / 0.4 rad/sec 
22.206 / 24.0 rad/sec 
4096 / 700 

20. 

90 

eta 

8 

0.0625 m 

3.0 inch 

408 kg/m? 

0 Pa in4 

180 sec 

Various 

AO 


Table 4.9: Series 7 Input Parameters 


Results 


Figures 4-22 - 4-25 compare the results between the FFT and SWS methods. The agreemeut 


between the methods is very good; the differences are on the order of the standard deviatious 


found in the previous section. 





Significant Wave Height 1.0 | 





Average Exposure 
Oo 9° 
o> ~ 
On O 
© oe) 






Tow Speed (kts 


oe So =r 


Figure 4-22: Comparison of SWS and FFT Methods: Average Exposure for 1.0 1m H, 
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Figure 4-23: Comparison of SWS and FFT Methods: Average Exposure for 3.0 4, 
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Figure 4-24: Comparison of SWS and FFT Methods: Threshold Met for 1.0 in H, 
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Figure 4-25: Comparison of SWS and FFT Methods: Threshold Met for 3.0 m H, 
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Chapter 5 


Future Efforts 


The Buovant Cable Array Autenna ATD is funded for FY 99 and FY 00. Should this proof- 
of-concept phase show that the floating antenna principles are feasible and cost effective. the 
Navy is expected to fund the development of a design which can stand up to the rigors of 
fleet deployment. Lincoln Laboratory is scheduled to test a prototype antenna in June 1999 
off the coast of Hawaii. Data will be collected on the electromagnetic performance of the 
altenha elements under varying conditions of source strength, source elevation angle. and source 
aziunuthal angle. In addition, the testing will document sea state and towing speeds so that 
performance can be related to the sea environment. 

Work at MIT will continue to support the floating antenna effort. Another graduate student, 
under the supervision of Professor Jerome Milgrain, will continue and extend the scope of the 
sunulation effort. Areas which will be investigated to improve the capability of the computer 


model include 


e 2 dimensional modeling of the sea surface 


e arbitrary angle between tow heading and direction of seas 


more detailed modeling of the hydrodynamics at the free surface 


higher order integration schemes 


validation against at-sea test data. 
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Appendix A 


User’s Guide to Antenna 


The User's Guide is broken down into two sections: 
e Input File 
e Output File. 


Section A.1 will go through the Antenna_ FFT input file line by line, and include program 
features or limitations associated with each paraimeter. At the end, features unique to the 
Antenna _ SWS input file will be addressed. Section A.3 will briefly describe the contents of each 
of the three output files, and recominended methods of viewing or extracting the information. 


The output files are the same for both prograin versions. 


A.l1 Antenna FFT Input File 


A.1.1 Header 


The input file for Antenna_ FF'T is shown in Appendix B. The first 7 lines are for informational 
purposes, and have no effect on the program. However, note that all variables are given in SI 
units. and that there must be 7 lines prior to the first variable defined. 


A.1.2 Run Identifier (run_ id) 


The run_id is a 7 character variable used to name the 3 output files. The 3 output files are 


named 


ce 





e ‘sn'// run_id 

e ‘ed’// run_id 

e °xl’// run_ id. 

(The // represents FORTRAN90 string concatenation.) Output files are located in the local 
directory from which the program is run. A sample pattern for run_id is run_id="001 FFT. 
Characters which ene would interpret as operators (/, +, -. *. etc.) should be avoided in 
run_ 2d if the Matlab© output option is selected. 


A.1.3 FFT Power of 2 (nnf) 


The FFT algorithin is designed to operate on arrays which have an integer power of 2 elements. 


nnf sets this power of 2. making 


N-sea = 2° (A.1) 


where 

N_sea is the number of sea frequency components and spatial locations. 

The nuinber of sea modeled in the spatial domain is generally much higher than the antenna 
length requires. however the speed advantage of the FFT method more than outweighs the 
penalty of carrying along the extra sea information. 


A.1.4 Number of Antenna Elements (num_of_ elem) 


This is the number of physical antenna receiver elements in the array. assuming that 


e the array is fully populated 


e the elements have zero space between them. 


A.1.5 Antenna Element Length (element length) 


This is the length of each antenna element. 


ffs 





A.1.6 Antenna Oversample Factor (M) 


Each antenna element is broken down into an integer nunber, M. of sections. These sectious 
correspond to the antenna spatial mesh size As. By enforcing the condition that As = Ax. the 


computer model constructs the spatial frequency meshing. Ab, such that 








] 
be = ihe 
ZENS ( ) 
fesea = As N sea (Aus 
i 
b = : 
. Pesea ee 


By invoking the dispersion relation. we convert from the 6 domain to the w domain. using 


Ww = /27gb (A.5) 


This process shows how the spatial meshing relates to the sea energy frequency domain 
meshing. As one can see, the fact that N_sea is constrained to be a power of 2 introduces some 


limitations on the mesh size As, since it is required that 
element_length = M As (A.6) 


A.1.7 Maximum Number of Samples (maz _num_of_ samples) 


This refers to the maximum number of “snapshots” taken for the Matlab© viewer, or the 
exposure data matrix. maz_num_of_ samples is used to dynamically allocate array storage for 


variables which are sampled periodically in time. 


A.1.8 Water Density (rho_ water) 


Water density in which the antenna operates. 


A.1.9 Antenna Density (rho_ antenna) 


The average antenna density, which is assumed to be constant over the eutire array. This is 


all approxination in that the prototype unit built by Lincoln Laboratory has shghtly higher 


ie 





deusity at the antenna element, and slightly lower density between eleinents. 


A.1.10 Antenna Diameter (diam_ antenna) 


This is the outer diameter of the hose-like structure encasiig the antenna. Element perforniance 
(vain) is controlled by whether this outer diameter is submerged; the diameter of the inner 
conducting slot anteuna affects the rotational invariance of the element, and the imaxinunu 


theoretical gain of a fully exposed element. 


A.1.11 Drogue Size (drog_ area) 


The drogue frontal area (in m*) determines the tension at the trailing end of the antenna. 
according to 


] 


as 5 Pwater 


TP U* drog-areaCD_n (A.D) 


A value of 0.01 m* was used for all simulations. 


A.1.12 Threshold Antenna Exposure (reqd_ ant _ exposure) 


As meutioned earlier. the antenna gain reniains steady until the element just becomes sub- 
inerged. In order to capture this performance. and introduce margin for capillary waves not 
modeled, the parameter reqd_ ant_ exposure has been defined. The reqd_ant_ exposure is the 
height above the water which must be maintained by the antenna to be considered exposed. 
The reqd_ ant_ exposure is given as a fraction of the antenna diameter. See Figure A-1 for a 
diagram representing an antenna dianieter of 1.0 and a reqd_ ant_ exposure of 0.04. (0.04 is the 


value used on all simulation results presented.) 


A.1.13 Tow Speed (U) 


The tow speed of the antenna is assuined constant with respect to time. Tow speed must be 
greater than zero (i order to generate tension); direction of tow is specified with the next vari- 


able. head_ seas. For the case of stern seas only, sea frequency components above omega-crit, 


omega_crit = - (A.8) 
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Figure A-1: Sample Threshold Antenna Exposure of 0.04 


are zeroed out to prevent antenna folding. If omega-crit is less than 1.5 rad/sec. then a warning 
is issued that excessive frequency components are being zeroed out. and program execution 
stops. 

A.1.14 Direction of Seas Option (head _ seas) 


The sea direction is limited to pure head seas or pure stern seas. head_ seas is a logical variable: 
.TRUE. indicates towing into head seas ( U < 0), and .FALSE. indicates towing with stern seas 
( U > 0). The sign of U is assigned internally by the program; U should always be positive 


when specified in the input file. 


A.1.15 Normal Drag Coefficient (CD_ 1) 


This represents the drag coefhcient for a cylinder in cross flow. A constant value of 1.0 has 


been assumed. 


A.1.16 Tangential Drag Coefficient (CD_ t) 


This is the skin friction associated with flow along the long axis of the antenna. A relatively 


smooth outer surface has been assumed in choosing this value, which has been assumed constant. 


Sl 





(0.0035 used in all simulations. ) 


A.1.17 Added Mass Coefficient (C_m) 


The added mass coefficient for a cylinder in cross flow has been set to 2.0. 


A.1.18 Simulation Time (duration) 


The total length of time simulated by the program. While many factors affect the actual tine 
necessary to perform the simulation (antemna length. number of mesh points. CPU speed. etc.). 
the ratio of actual time to simulation time is roughly 5:1. (For the SWS version. the ratio is 


roughly 50:1.) 


A.1.19 Time Margin (margin) 


The tine step size used by the program, t_ incr, 1s given by 


. _ As 
tincr = margin— (A.9) 
C 
where c is the string wave propagation velocity 
i 
Pant "7 


ald 

Tis the string tension 

Pant 18 the antenna density 

D is the antenna diameter. 

The tension used initially to estimate the time step is the calm water (local fluid particle 
velocity. u = 0) tension at the tow point. When dynamic_ tension = .TRUE., then the actual 
tension varies with x and ¢t as the actual wetted surface area changes, and as the local fluid 
particle velocity. u. evolves along the length of the antenna. A sample plot of inaxinum tow 
point tension vs. significant wave height using the parameters of Table A.1 is provided in Figure 


A-2. As can be seen in the figure, the maximum tow point tension in a 3.0 m significant wave 
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Parameter Description 


Method SWS 

ye various 
Devel 1.0 

Wmin 0.4 rad/sec 
Cee 24.0 rad/sec 
Deen 700 

L 29 1 
Nelem 30 

eres. Oe sara 

nu 8 

As 0.0625 m 
Diameter a. Oinen 
pe 408 kg/m 
ee 0 Pa m7 
Duration 180 sec 

U Duns 
Netnresh 40) 


Table A.1: Tension hiput Paraiueters 


200.0 
180.0 
160.0 4 
140.0 
120.0 
100.0 

80.0 


60.0 5 


Maximum Tow Tension (N) 


40.0 4 


20.0 





0.0 , 
0.0 0.5 1.0 es 2.0 Ze oe a5 
Significant Wave Height (m) 


Figure A-2: Variation of Maximuin Tow Point Tension vs. Significant Wave Height 
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height sea is about triple that of the calin water value. This would lead to a required value for 
margin of ee or about 0.58. In practice, a value of 0.25 to 0.33 was used: however. no detailed 
effort was made to determine the maximum value allowed for margin. 

For the case of dynamic_ tension = .FALSE., the calm water tow tension is also the max- 
imum tow tension. (Tension is determined entirely by a relative fluid velocity of U. aud the 
wetted area corresponding to the cali water submergence depth.) Thus a value much closer to 
1.0 is possible for margin. 

When the antenna is modeled as having lateral stiffness G.e. fler_ stiff_ coef > 0.0). the 
variable margi is used to empirically determine a satisfactory time increment. Time step size is 
quite sensitive to aitenna stiffiess, aud margin must reduced an order of magnitude for typical 


designs considered (e.g. fler_ stiff. coef =100 Pa in‘). 


A.1.20 Snapshot Output Interval (pic_ time) 


This variable allows the user to control the time interval between snapshots. which are recorded 
to file ’sn’//run_ id. These snapshots record both the sea and antenna elevation data along the 
entire length of the antenna, at a resolution of As. A Matlab© viewer script file is included 
in Appendix D, which allows the user to step through the snapshots frame by frame. pic_ time 


has been set to 1.0 second for all run results presented. 


A.1.21 Time of First Snapshot (t_ first pic) 


The time to begin recording snapshots is user defined. This is useful to ignore transient behavior 
at the start of the simulation, or to save only a representative sample at the end of a lengthy 


sluiulation. 


A.1.22 Statistics Output Interval (stat time) 


This variable allows the user to control the time interval between exposure data matrix records. 
which are recorded to file ’ed’//run_ id. While stat_ time is independent of pic_ time, a value of 
1.0 secoud was also used for this variable for all simulations. Selection of a stat_ time excessively 
siiall results i subsequent records containing essentially the same information, and contributes 


little to the statistics of the sainple. 
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A.1.23. Threshold Number of Elements (num_ of elem_ thresh) 


This variable allows the user to specify the lowest number of exposed antenna elements necessary 
to achieve the desired design antenna array gain. At each stat_ tame. the number of exposed 
elements is compared to num_of_elem_ thresh; the result (a “1” when the threshold is met. 
“0” otherwise) is written at the eud of the summary output file. “xl’//run_ id. along with anv 
other variables tracked over the course of the simulation. Currently. the variables tracked over 


tlie are 
e tow point tension (tow_ tension( )) 


nuuber of elements exposed (is exposed( , )) 


threshold met (thresh met: “O” or “1”, as above) 


e counterdetectiou leugth (length_ totally exposed( )). 


A.1.24 Tension Model Option(dynamic_ tension) 


(See discussion under Section A.1.19.) 


A.1.25 Snapshot Output Option(matlab output) 


The snapshot output file (sn’//run_ id) is optional, since it does not contain any unique sta- 
tistical information, and requires a relatively large amount of disk space. (Average file size is 
about 1.4 MB for a 180 second simulation.) The primary benefit of examining the snapshot 
output file is to learn the physical behavior of the antenna in a seaway. In addition. the output 
is useful to detect when the time step size is too small, which is characterized by rapid un- 
physical antenia oscillations. matlab_ output=.TRUE. means that suapshot output. is enabled, 


FALSE. means snapshot. output is disabled. 


A.1.26 Significant Wave Height (sig wave_ht) 


This variable sets the sea’s significant wave height, a parameter used to describe the sea sever- 
ity in both the Bretschneider and Ochi spectral models. Appendix E correlates sea state to 


significant wave height for a fully developed sea. 
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A.1.27 Development (development) 

As discussed in Section 2.1.2, this variable allows the user to shift the sea spectral energy modal 
frequency up (for developing seas) or down (for decaving seas). This parameter applies only to 
the Bretschneider spectral model. 


A.1.28 Spectral Model (spectral_ model) 


Two spectral model subroutines are available: the two parameter Bretschneider spectrum. and 
the Ochi “most probable” six parameter spectrum. (“1” means Bretschneider. “2” means Ochi.) 


The Bretschneider spectral model was used for all simulations presented. 


A.1.29 Antenna Lateral Stiffness (flex stiff_ coef) 


This variable represeuts the resistance of the antenna to bending, and is given by 
iilewostt}| f_coey = HT (A.11) 


where 
Fis the modulus of elasticity 


T is the area inonmient of imertia. 


A.2 Input File Features Unique to Antenna _ SWS 


The differences in the Antenna _ SWS input file pertain to the discretization of the sea circular 
frequency components. 
A.2.1. Number of Sea Frequency Components (num_ of_ freqs) 


This variable is the number of discrete sea frequency components in the spectrum. num_ of_ fregs 
= 700 was used for the simulations, however as few as about 100 components provide adequate 


sea representation. The frequency components are evenly spaced in the w domain. 
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A.2.2 Frequency Coverage (omega_ min, omega_ maz) 
The frequency components w, are eveuly spaced fro omega_ mun to omega_maz. For the 


simulations presented omega_ min = 0.4 rad/sec and omega_ maz = 24.0 rad/sec. 


A.2.3. Spectral Model (spectral model) 


One extra option is added: for spectral_ model = 3, the sea surface is described by a pure 


harmonic wave of amplitude s7g_ wave_ ht and circular frequency omega_ maz. 


A.3 Output Files 


A.3.1 End-State Summary File (’xI’//run_ id) 


This file provides 


e a printout of forces, velocities, and positions at the end of the run 
e a record of the input parameters used for the simulation 
e averaged statistics results, and 


e a time history of variables tracked by the program. 


®) 


This file is suitable for viewing with a spreadsheet application, such as Microsoft Excel~ . 


A.3.2. Exposure Data Matrix File (’ed’//run_ 1d) 


This file, whose data structure is given in Table 4.1, contains the raw exposure data information 
for each anteuna element at each tine sampled. This allows one to analyze statistics across 


tie or space. This file is easily analyzed by a spreadsheet application. 


A.3.3 Snapshot File (’sn’//run_ id) 
®) 


The snapshot file can be viewed with a user written code, or by Snaps.m, a simple Matlab 
routine provided in Appendix D. It provides views of the antenna as it floats on the surface at 


two scales: one view is large enough to see the most extreme waves, one zooms in on the sea 


7 





surface by using the local water elevation as the zero for the z axis. The latter view also shows 


the antenna ceutral axis. the antenna upper limit. aud the antenna lower liunit. 
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Appendix B 


Input File for Antenna FFT 


Input File for Antenna_FFT Simulation Program 


SI units throughout; useful conversions follow 


1 inch 
ie Tooc 


1 knot 


OCOLFFT 
13 

25 

0.5 

8 

300 
1025.9 
408.0 
0.0762 
2001 
0.04 
2.9722 
- TRUE. 
1.0 
0.0035 


oO 
© 


0.0254 m 
0.3048 m 
0.5144 m/s 


run_id 

nnf 

num_of_elem 
element_length 

M 
max_num_of_samples 
rho_water 
rho_antenna 
diam_antenna 
drog_area 
reqd_ant_exposure 
U 

head_seas 

CDin 

CD_t 

C_m 

duration 

margin 

pic_time 

t_first pre 


stat_time 


7 characters long 

# of components. N_sea=2*¥*nnf 

# of actual slotted antenna elements 
physical slotted antenna length (m). 
# of mesh points/element_length. 

to dimension time arrays 

seawater at 15 deg C (kg/m°3) 

density of antenna (kg/m73) 

antenna outer diameter (m) 

drogue area normal to flow (m2) 
threshold vertical exposure for reception. (x/diameter) 
tow speed (use POSITIVE #) (m/s) 

TRUE: antenna is towed into head seas 

normal flow drag coefficient 

fluid tangential drag coefficient (p. 50 Berteaux) 
added mass coefficient 

total simulated run time (sec) 

saftey factor applied to time step size 

snapshot interval for Matlab viewer (sec) 

time of first snapshot (sec) 


statistics interval (sec) 
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20 


TRUE. 
- TRUE. 


20 
ao 


num_of_elem_thresh 
dynamic_tension 
matlab_output 
Ssig_wave_ht 
development 
spectral_model 


Plexestattecoet 


# of exposed elements reqd for specified gain 

FALSE: steady-state calm water model to find tension 
TRUE: Matlab viewer output is ON 

significant wave height (m) 

modal frequency shift factor 

i=Bretschneiaer; 2-0chi 


E*I] (Pa m4) 
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Appendix C 


Input File for Antenna _ SWS 


Input File for Antenna_SWS Simulation Program 


SI units throughout; useful conversions follow 


Pinch 
1 foot 


1 knot 


QO00_SWS 
25 

0.5 

8 

100 
300 
1025:.9 
408.0 
0.0762 
OF01 
0.04 
2.0022 
, PRUE. 


.0035 


— 


0.0254 m 
0.3048 m 
0.5144 m/s 


run_id 
num_of_elem 
element_length 

M 

num_of _freqs 
max_num_of_samples 
rho_water 
rho_antenna 
diam_antenna 
drog_area 
reqd_ant_exposure 
U 

head_seas 

CD_n 

CD_t 

C_m 

duration 

margin 

pic_time 

CoLirste pic 


stat_time 


7 characters long 

# of actual slotted antenna elements 
physical slotted antenna length (m). 
# of mesh points/element_length. 

# of discrete freq components 

to dimension time arrays 

seawater at 15 deg C (kg/m°3) 
density of antenna (kg/m73) 

antenna outer diameter (m) 

drogue area normal to flow (m°2) 
threshold vertical exposure for reception. (x/diameter) 
tow speed (use POSITIVE #) (m/s) 

TRUE: antenna is towed into head seas 

normal flow drag coefficient 

fluid tangential drag coefficient (p. 50 Berteaux) 
added mass coefficient 

total simulated run time (sec) 

saftey factor applied to time step size 

snapshot interval for Matlab viewer (sec) 

time of first snapshot (sec) 


statistics interval (sec) 


ul 





20 

~ RUE. 
. TRUE. 
0.40 
15.0 
2.0 
1-9 


num_of_elem_thresh 
dynamic_tension 
matlab_output 
omega_min 
omega_max 
Slig_wave_ht 
development 
spectral_model 


flex_stiff_coef 


# of exposed elements reqd for specified gain 

FALSE: steady-state calm water model to find tension 
TRUE: Matlab viewer output is ON 

lower cutoff freq of sea spectrum (rad/sec) 

upper cutoff freq of sea spectrum (rad/sec) 
Significant wave height (m) 

modal frequency shift factor 

1=Bretschneider; 2=Ochi; 3=pure sine wave 


E*I (Pa m4) 
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Appendix D 


Source Code Listing for Snaps.m 


' snaps displays antenna and sea surface as a series of snapshots 
! displays series of antenna snapshots 


' N=number of antenna elements 


' M=number of snapshots taken 


' number of rows in file M+1 


' number of cols in file = 2*N+1 


' File format (ASCII) 


SGOL ol: 

{ zero (a dummy), time 1, ... , time M 

micol 2: 

| mole zil Pate tame t. 2220, 2(1) at time M 
P coL Nai: 

' mite Zin) at time 1. ... , Z(N) at time M 
' col) N+2: 

| x(1), h(€1) at time 1, ... , h(N) at time M 
! col 2*N+1: 

ieee xo n(N) at time 1, ..., HCN) at time M 
clear 


load snOOO_SWS -ascii 
(M,NJ]=size (sn000_SWS); 


N=(N-1)/2; 
M=M-1; 
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t=sn000_SWS(2:M+1,1); 
x=sn000_SWS(1,2:N+1); 
z=sn000_SWS(2:M+1,2:N+1); 
h=snO00_SWS (2:M+1,N+2:2*N+1) ; 
diam=sn000_SWS(1,1); 
y=z-h; 
for wet. N 
zz(j)=0.0; 
thresh1(j)=0.04*diam; 
thresh2(j)=-0.04*diam; 


end 


' find max/min z for consistent plot scales 
Z_min=min(h); ' z_min 1S a row vector 
Z_m1n=min(z_min) ; ' z_min 1S now a scalar 
z_max=max(h); 
Z_max=max(z_max); 
figure(1); 
for j=1:M 

ej! 

pause; 

end 

subplot (2,1,1); 

pion Gc (1,2), x,nC)52), *,22); 

title([’Time’, blanks(3),num2str(t(j))J); 

ylabel(’z, h (m)’); 

axis((0.0 x(N) z_min z_max]}); 

subplot (2,1,2); 

plcuGc wit). 6: ’,x, thresh! ,’b-" x, thresh2?,’b-’,x,zz,’e-") : 

axis([0.0 x(N) -.05 .05)); 

xlabel(’x (meters)’); 


ylabel(’y =z-h (m)’); 


end 
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Appendix E 


Sea State Data for Fully Arisen Seas 





Significant Modal Wind Sea 
Wave Height Frequency Speed State 


(an) (rad/sec) (knots) 
0.1 3.96 4.2 
0.2 2.80 9.9 
0.3 2:29 3 l 
0.4 1.98 8.4 
Ore eee 9.4 
0.6 eG? 10.3 
One 1.50 Piel 
0.8 L.A0 11.9 2 
0.9 ioe 226 
1.0 iz 13.3 
el ie20} 13.9 
ee ore 14.6 
is 1.10 2 o 
ied 1.06 oe 
eo 1.02 6:3 
Leb 0.99 16.8 
7 0.96 eo 
1.8 0.93 as 
1.9 0.91 oc 
220 0.89 18.8 4 
D2 0.85 19.7 
hy 0.81 20.6 
2.6 0.78 24 
2.8 0.75 Dee 
3.0) Oee2 23.0 o 
BD Ce 23.8 
3.4 0.68 2400 
3.6 0.66 Do 6 
Be 0.64 Pao 
4.0 0.63 26.6 
Ans 0.59 28.2 
5.0 0.56 29.7 
ED Ors 31.2 
6.0 0.51 B20 
So 0.49 Sore T 
eo OAa7 Oo 
C5 0.46 36.4 
8.0 0.44 Deo 
8.5 0.43 38.8 
9.0 0.42 39.9 


Table E.1: Sea State Data for a Fully Arisen Sea 
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